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Abstract: We study the action for a non-Abelian charged particle in a non-Abelian 
background field in the worldline formalism, described by real bosonic variables, leading 
to the well known equations given by Wong. The isospin parts in the action can be viewed 
as the Lagrange multiplier term corresponding to a non-holonomic constraint restricting 
the isospins to be parallel transported. The path integration is performed over the isospin 
variables and as a result, the worldlines turn out to be constrained by the classical solutions 
for the isospins. 

We derive a wave equation from the path integral, constructed as the constrained Hamil- 
tonian operator acting on the wave function. The operator ordering corresponding to the 
quantum Hamiltonian is found and verified by the inverse Weyl transformation. 
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1. Introduction 

A starting point in the usual way of quantization for a system is to adopt its classical 
system as a fundamental framework. Interestingly, different from the case of an electrically 
charged particle, the dynamics of a particle with non-Abelian charge was quantized directly 
from analogy and extension of its quantum theory since a non-Abelian charged particle had 
not been observed in a classical context. If a classical dynamics of a particle is considered 
as the most fundamental basis to construct a quantum theory, to find a classical picture of 
a particle with non-Abelian charge can be not only an interesting quest but also a crucial 
part for understanding its quantum nature. 

A set of classical equations of motion for a non-Abelian charged particle was first 
derived by Wong from the quantum action in which a non-Abelian charged particle is 
described by a Dirac field [1]. These are called Wong's equations and include the parallel 
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transport equation describing an isospin in addition to a non-Abelian extension of the 
Lorentz force equation. 

The purpose of this paper is to find a classical picture for a non-Abelian charged 
particle in a background field. We start from a classical action for a non-Abelian charged 
particle and construct the path integral. Many actions producing the Wong equations 
have been proposed [2-5]. Since we are interested in a purely classical picture, we write 
the action using only variables valued in the real numbers. In particular, all our c-numbers 
commute with one another, so upon quantization they will become bosonic variables. In [2], 
using the constrained Hamiltonian formalism, a wave equation governed by the constrained 
Hamiltonian operator was written down though the operator orderings were unspecified. 
We find that the wave equation gets an additional term which breaks gauge invariance. 
This will be verified by an inverse Weyl transformation. The derivation brings forth certain 
issues of operator ordering. 

As is usually done, for instance in [2, 3], we first consider the equations given by Wong 
as a starting point to eventually obtain a quantum description for a non-Abelian charged 
particle. The first step is to obtain the classical action producing these equations. As 
mentioned above, we use only commuting real variables for all degrees of freedom. Next we 
quantize this action using the constrained Hamiltonian formalism due to Dirac [6, 7]. Once 
the Hamiltonian is obtained, we construct the path integral in the worldline formalism, 
and derive the constrained wave equation. This equation is compared with the one derived 
in canonical quantization earlier, discussing an operator ordering issue [8]. 



2. Classical action for a non-Abelian charged particle 

A classical point particle is characterized by its position in spacetime, and if it has 
non-Abelian charge, also by a vector in some internal space, corresponding to some repre- 
sentation of the gauge group. For now, we will take this vector I a to be in the Lie algebra 
of the gauge group, and refer to it as the 'isospin' of the particle. We can also take the 
internal vector to be in the fundamental vector space of the gauge group, we will discuss 
that later. 

Then the following well-known pair of equations describes the dynamics of a classical 
non-Abelian charged particle, first derived from the Dirac equation by Wong [1], 



d 

m— 

dt 



i a + gf abc A b ^ = 0, (2.1) 
+ gI a FlX = 0. (2.2) 



(_^)l/2j » ^ 



The equation of motion for the isospin I a {t) is a parallel transport equation along the 
trajectory of a particle, x^(t), with the connection T acfl = gfa^A 1 ^. The gauge group will 
be taken to be compact semi-simple, and the structure constants j a \, c will be chosen as real 
and totally antisymmetric. The equation of motion for x^(t) is a non-Abelian generalization 
of the Lorentz force equation. The parameter t is a parameter along the worldline of the 
particle, and a dot denotes differentiation with respect to t. Let us consider the following 
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action, 



J dt (-mi-xHtfl 2 + J a {i a + gfabcA^r 



(2.3) 



The variables appearing in this action take only real numbers as values, and will thus 
correspond to real bosonic variables after quantization. The background metric is r]^ = 
diag(— 1, 1, • • • ,1). We have not included a term for the dynamics of the gauge field, as we 
ignore the backreaction of the charge on the field. The Euler-Lagrange equation obtained 
by varying J a (t) is the equation of motion for I a (t), 



Q= ddL L _dI L = _ ia _ , 
dtdJ a dJ a yJ 



(2.4) 



Similarly, varying I a (t) gives a parallel transport equation for J a (t), 





— = J a + gfabcA^X >*J C . 



dt dl a dl a 



(2.5) 



The equation of motion for x^(t) is a non-Abelian version of the Lorentz force equation as 
may be expected, 



m 



d_ 

dt 



(-x v x v yi 2 



+ gK a F^x» = 0, 



(2.6) 



where K a = f a bcJ b I c = (Jx I) a . It is easy to check that K a also satisfies the same parallel 
transport equation as for J a and I a . Therefore, the action Eq. ( |2.3[ ) leads to the same set 
of equations Eq. (^) given by Wong. 



K a + gfabcAix»K c = 0, 



m- 



dt 



{-x v ± v fl 2 



+ gK a F^x» 



0. 



(2.7) 
(2.8) 



We note that the action given here is essentially the same as the one given by Balachandran 
et al. [2], but written in terms of real variables I and J rather than the complex variable 
a = I a + iJ a . 



2.1 Gauge invariance 

The action is invariant under the following set of infinitesimal transformations, 

s[i] a = 9fabcA b I C , 
d[A}- =gf abc A b A^-d II A a , 



S[FL] = gfabcA'K 



(2.9) 



where A a are infinitesimal. The same transformation rule as for I a applies to J a and 
K a . Under these transformations the parallel transport equation Eq. ( |2.4[) transforms 
covariantly while the generalized Lorentz force equation Eq. ( |2.6| ) is invariant. For the 
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parallel transport equation, 

5[i a +gf abc {A b ^I c )] = + gf abc 5[A^I c + gf abc Alx»5[I] c 

= gfakiA k i l + gfakiA k i l + gfa bc (gfbkiA k A 1 ^ - A b )r 

+gfabcAl^{gf ckl A k I l ) 
= gfakiA k (i l + gfi bc A b l x»I c ), (2.10) 

where we have used the Jacobi identity. Similarly, Eq. (|2.6| ) is invariant under the gauge 
transformation since F£ v and K a transform covariantly. 

The action has a further invariance under reparametrizations of the worldline. If t is 
replaced by some smooth function r(t) with f 7^ , we must make the replacements 

d d 1 d ,„ _ _ s 

dt^dr = Tdt, — ^ — = — -. (2.11) 
at dr f dt 

It is easy to see that the action Q2.3] ) remains invariant. 



3. Constrained Hamiltonian formalism 



The Hamiltonian plays a central role in the time evolution of physical quantities. However, 
the presence of gauge symmetries makes this a constrained system — the variables in 
the Lagrangian are not all independent. These symmetries include invariance under the 



transformations of Eq. (2J3), as well as invariance under reparametrizations of the worldline. 
We will follow the treatment due to Dirac [6, 7], which we describe briefly here. 

The Lagrangian is a function of (q,q). When the Hamiltonian is obtained by the 
Legendre transformation, a variable q is converted to the canonical momentum p defined 
by P = §r- This definition may lead to some constraints, which are relations among the 
canonical variables, of the form 4>m(q,p) = 0. These are added to the original Hamiltonian 
Hq = pq — L with Lagrange multipliers u m (q,p), giving a new Hamiltonian, 

H(q,p) = H (q,p) + u m {q,p)4> m (q,p). (3.1) 

Next we require that the constraints should not change with time (the consistency con- 
ditions). These conditions can lead to new constraints as well as relations among the 
Lagrange multipliers. 

O = m = [<A m ,#o]+w™[<Am,M (3.2) 

where the square brackets signify Poisson brackets. That is, these equations may give a 
new constraint (a secondary constraint) or specify unknown coefficients u n . 

(a) If a new constraint is found, it is again applied to the consistency condition in 
order to see whether it gives a new further constraint or restricts the coefficients u n . This 
algorithm ends if no more new constraint comes out. 

(b) If relations are found among the Lagrange multipliers, u n , general solutions for 
u n to Eq. ( |3.2[) are written in the form u n = U n + V n , where U n are particular solutions 
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and hence specified, and V n are homogeneous solutions, which are in a form of linear 
combination of the solutions with arbitrary coefficients v a , V n = v a V™- 

Any constraint which has vanishing Poisson brackets with all constraints is called a 
first class constraint and all other constraints are called second class constraints. The 
consistency conditions ensure that the Hamiltonian has vanishing Poisson brackets with 
all the constraints. Ultimately all the constraints including second class constraints need 
to be imposed on the system in order to keep only independent variables in the theory. 

In order to do this, we first define Dirac brackets using all the second class constraints, 
denoted by Xa, 

[F,G]* = [F, G] - [F, X a]C af3 [x0,G], (3.3) 

where C is the matrix given by the matrix elements C a p = [x a , Xp] an d the inverse elements 
C°P = [C~ 1 ] a p. The Dirac bracket, like the Poisson bracket, is antisymmetric and satisfies 
the Jacobi identity, but the Dirac bracket of any quantity with a second class constraint is 
zero, 

[ Xa ,F}* = 0. (3.4) 

This implies that by replacing the Poisson brackets by Dirac brackets [F, G]*, we can set 
all constraints, including the second class constraints Xa, equal to zero at the end of a 
calculation. 

3.1 Worldline formalism 

For technical convenience we will use the worldline formalism to write the action [9, 10], 
and then to compute the path integral, for the particle with non-Abelian charge. The 
action in Eq. ( |2.3D can be rewritten in the worldline formalism by introducing a Lagrange 
multiplier h, 

S = y^(L + L/), (3.5) 

where 

L = \h~ l x 2 - h^-, Lj = J a [i a + g(A fl x fM x /)"]. (3.6) 



The introduction of h makes the action in the worldline formalism invariant under the 
reparametrizations 

t -►/(*), h^h/(df/dt). (3.7) 
The equation of motion for h is 

x 2 + h 2 m 2 = 0, (3.8) 
while the other equations of motion are now 



i a + gfab C A b u r±» = o, (3.9) 



d 



-[h-H^+gK^F^ = 0, (3.10) 



dt 



where K a = f abc J b I c = (J x I) a as before. Using Eq. Q in Eq. (|3~To|) , we recover 



Eq. (2.2). Alternatively, if Eq. ( |3.8[ ) is plugged back into the action, the previous Lagrangian 
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Eq. (|2.3j ) is recovered. The canonical momenta can be calculated from Eq. ( |3.6D , 



Pj = 0, (3.11) 

Pf = J a , (3.12) 

P h = 0, (3.13) 

P» = h-H^ - gA°K a , (3.14) 

from which we can read off the constraints immediately, 

<Pi=Pj^ o, 

0" = Pf - J a RJ 0, 

0h = Ph ~ 0, (3.15) 

where the symbol ~ indicated weak equality, i.e., equality on the submanifold defined by 
constraints. Now we write down the canonical Hamiltonian Hq, 

H = P^ + Pfi a + Pjj a + P h h - (L + Lj) 

= ±[(P^ + gAlK a f + m 2 \ + PU a + {Pi-J a )i a + Phk (3.16) 

where it has been assumed that the time derivatives of the coordinates, e.g. I, are expressed 
in terms of the canonical variables, in principle. If the constraints Eq. ( |3.15|) are added to 
Hq, the new Hamiltonian H\ remains in the same form as before, 



H x = H + C a x P a j + c%{pf - J a ) + C h P } 



h 



-I 

2 [ 

+P h [C h + h(p,g)]. (3.17) 



-[(P, + gA^f + m 2 ] + P?[C? + J a (p, q)} + (Pf - J a )[Cf + i a (p, q)} 



The consistency conditions require that the constraints have vanishing Poisson brackets 
with Hi, so writing IT M = (P^ + gA®K a ) , we find 

= 0? = [#,#!] 

= [Cf + r (p, q)\ + 2 5 (/ afec 4/ c )(/ l /2)m t , (3.18) 
O = 0| = [<f>%,Hi] 

T a (r> nW -On(f, ^ , 



-[C? + J a (p,g)] - 2 5 (/ afec ^J c )(/ l /2)lT t , (3.19) 



= 0,,= [i\,fTi] 

= -i[(P M + 5 ^) 2 + m 2 ]. (3.20) 

The first two relations allow us to eliminate Cr? + I a (p,q) and Cf + J a (p,q) in i?^. The 
last equation is a secondary constraint, which we denote by 03 , 

03 = (P„ + 5^^ a ) 2 + m 2 « 0. (3.21) 
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Now we can rewrite H\ with the unknown coefficients eliminated, 



-(U^ + m 2 ) + (C h + h)P h 



+ 'i:[-2g{f abc A b li J c )WP a j - 2g(f abc A b fl I c W(P? - J a )\ 



(3.22) 



where we have defined 11^ = P^ + gA a ^K a . We can see that 4>3 identically commutes with 
Hi without giving any further constraint or condition for the coefficients, 

4>3 = [(/>3,Hi] 



h 



2 3 / flk ^J c n"[n 2 + m 2 , P a j] - 2gf abc A b J c W"[U 2 + m 2 , Pf - J a ] 



Agf abc A^J c U^[U u , P a j] - 4gf abc A^I c U^[U u , Pf) 
Also, by the following consistency condition, must vanish, 

h = [h, Hi] = (C h + h)[h, P h ] = C h + h. 
Thus the final expression for the Hamiltonian is 



0. 



(3.23) 



(3.24) 



H = h 
with the constraints 



-(IPII^ + m') - gf abc A^J c P] - gf abc A^P{Pf - J a ) 



pa 



4>3 
01, 



I = pf - r « o, 

a i^a\2 



{P^ + gA^K 
Ph-0. 



m 



+ hP h , (3.25) 



(3.26) 
(3.27) 
(3.28) 
(3.29) 



and 



can be recombined into one first class and two second class con- 



Note that 

straints. The term within the square brackets in Eq. (3.25) is the first class linear com- 
bination of these constraints. We may call this term <pQ. Then we can choose <f)f and ^ 
as the independent second class constraints, and 4>q and 4>h as the independent first class 
constraints. In other words, the constraint in Eq. (3.28) is replaced by <f>o ~ 0. We note 
that the Hamiltonian of Eq. ( |3.25j ) is also a first class constraint by construction. If the 
gauge freedom in terms of arbitrary coefficients is considered, H does not explicitly show 
the gauge freedom for the first class secondary constraint <f)Q. For this, one can use the 
extended Hamiltonian, where the number of unknown coefficients explicitly matches with 
that of the first class constraints, 



H E = H + 



(3.30) 



Recognizing h is arbitrary we see that the number of the arbitrary coefficients is equal to 
that of the first class constraints. 
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We can calculate the Dirac brackets using the matrix 

C12 = ^ah = — Cfl) (3.31) 

and set the second class constraints to zero due to the property [xa>-^1* = where Xa is 
the second class constraint and F is an arbitrary quantity. The constraint Ph ~ 0, along 
with the gauge fixing condition h = A, which we now set, allows us to remove a pair of 
variables (h,Ph) [9]. Then we find 

H E = C(U^ + m 2 ), (3.32) 

as also given in [2]. One can see whether we use H or He the path integral will be in the 
same form as long as all the corresponding gauge conditions to the first constraints are 
applied for it. 

4. Path integral 

A general form of the path integral for a constrained system is given in [7, 11]. 

* A II 5( X *)(det{ Xa , xp]) l/2 II WaWGa) n det[G a , <f> b ]ei s ^ , (4.1) 



1Z 

t,a t,a 



where 4> a and Xa are a first and a second class constraints, respectively and G a is a gauge 
fixing condition for <j) a . Here 5(x a ) is responsible for the second class constraints and 
(det^a, X/?]) 1//2 f° r their Jacobian factor. After suitable gauge conditions for the first class 
constraints are chosen, a gauge condition and its corresponding first class constraint form 
a pair of second class constraints. Thus the delta functions corresponding to the first 
class constraint and the gauge fixing condition, 5((j) a ) and 5(G a ), respectively, are similarly 
incorporated together with the Jacobian factor det[G a , <pb\ in the path integral. The gauge 
fixing constraints need to be introduced in order to fix the arbitrary coefficients on the 
corresponding first class constraints in the Hamiltonian. Applying the formula Eq. ( |4.1| ) 
to our case, we find that the path integral becomes 

J™ d\ J VI P^/^P^Mlli - ^^,XoM^]e ilW '" lp " 1 '^ P? . (4.2) 

The Jacobian factor (det[x«, Xp]) 1 ^ 2 = 1 an d the delta functions for the second class 
constraints have been used by setting Pj to zero and replacing J a by Pf . The integration 
over the gauge choice A for the first class constraint Ph « implies that the gauge freedom, 
i.e., a reparametrization invariance, is factored out [9]. «S^ n is 

%in = Pvfl? + Pfi a + PjJ a + Phh - H 

= IL^ - gA^K 11 + Pfi a - hcp . (4.3) 

We already noted that the Hamiltonian itself is a constraint and hence zero. (j)o can a lso 
be placed in the same footing. One can express the delta function for c/>o in the following 
form 

«5[0o] = — / ^uje 1 ^^ . (4.4) 

27T ./ 
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Hence, the delta function for cpo can be combined into the Hamiltonian by redefining 
h — to — > h, 

~ J QheU^W* (4.5) 

Therefore, including the delta function for the constraint (f>o is redundant. 

After the redefinition of h we put the gauge condition S[h— A]. Then we notice that Jz^ n 
becomes different from the original Lagrangian Lq + Lj, Eq. fl3,6|) , only by — ^ (II — \x\ . 

%in = + Pfl a + P a jj a + P h h-H 

= - gA a ^K a + Pfi a - |n 2 - ^m 2 

= _ K n _ W + ^ _ ^ m2 + ja[ia +a{A ^ x /r] 

+{Pf - J a )i a + P]j a . (4.6) 

where x^ = Xn+ ^ t Xn ,etc, so Jz?j n has been expressed only by canonical variables (p, q), i.e., 
Jzfj n = J£i n {jp,q). After the terms (Pf — J a )I a and Pj J a are removed by the second class 
constraints with the Dirac bracket formalism, the Lagrangian becomes as follows. 

Se^—U- ±i) 2 + ^x 2 - \m 2 + Pf[i a + g(A^ x /)"]. (4.7) 

According to the definition of the momentum P„, the first term would be zero in continuum 
limit. However, one cannot simply set it equal to zero since x^ = Xn+ ^ t Xn is treated as 
independent variables in the path integral. This implies that the path integral directly 
expressed from the given Lagrangian may be different from one reconstructed from the 
Hamiltonian. Our interest in this paper is in the integration over the internal variables 
(Pf, I a )- In the next subsection we see that a certain gauge choice can allow us to integrate 
over the internal degrees of freedom. 

4.1 Gauge condition for 0o 

Here we choose a gauge fixing condition \o f° r <fio- The gauge condition comes into the 
path integral as the delta function 8[xo] with the Jacobian [4>q, xo\- We are interested in the 
case that the gauge fixing does not spoil form of the isospin part so that later we are able 
to integrate over the internal variables separately. We choose the gauge fixing condition, 

Xo = x° - t = 0. (4.8) 
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It is easy to check the Poisson bracket [<fo>Xo] is non-zero 

1 



)>Xo] 



-(IT + m 2 ) + Cf + C^, x v - t 



-{U 2 + m 2 ),x° -t 



1 



= ^x2U^[P II + gAf l K a ,x } 



2 

= -n°. 



(4.9) 



By a change of variables Il^—^x^ — >■ P M , the Jacobian factor [0q, Xo] becomes (— P° — j^x°) . 
The relevant path integral can then be written as 



/ 



-P° - -x° 



1 



e h ■ 



(4.10) 



9h8[h - \]@x°6[x° - ^^^ei/^l-l^+^^-l^+^^+s^^x/)"]} 

;/rft[-|(Po) 2 + ^(io) 2 ] 

„-J ^/ M y[/ l _A]^x i e^*U*J-V+i?^+9(^x7)-]} is0=tj 

where the Gaussian integrals over P M have been taken. 
4.2 Integration over internal variables 

In the last subsection we have chosen the gauge fixing condition for 4>q. Since this gauge 
fixing does not change structure of the isospin part, we are able to compute an integration 
over isospin variables. Taking the relevant part of the action in a time-sliced form, 

N+l 



Si = ^2 Pl,n[In ~ In-1 + 9fab C A^ n (x n — X n _ 1 )i n ] 
n=l 
N+l 

= J2 P UK C In-I C n-l), 



(4.11) 



n=l 



where M® c = 5 ac + g f a b c A^ n (xn — Xn-i)- The integral over Pf n in the path integral 
produces a delta function, 



J dP^e^J^-i^) = 27 rM(il« - 



(4.12) 



What finally remains in the path integral is a product of all the delta functions and after 
integration over I a it reduces to a single delta function, 

{2ttK) n+1 j dI a N --- dI?5(M N+l I N+1 - I N )S(M N I N - I N _i) • • • 6(M 1 I 1 - I ) 

= (2irh) N+1 J dI a N _i • • • dI15{M N M N+l I N+l - /jv-i) • • • 6 (Mill - I ) 

= (2irh) N+1 5(Mi ■ ■ ■ M N M N+ iI N+ i - I ). (4.13) 
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The matrix product M\ ■ ■ ■ MjyMjy+i is a finite parallel transport operator backward in 
time, but it can be placed in time order. Let us define (A n ) ac = gf a b c A^ n (xn — x^_ 1 )/At n , 
where At n = t n — t n -\. Then using the antisymmetry of A n , we can write 



[Mi • • • M N M N+1 ] ab I b N+1 

[(1 + AiAii) • • ■ (1 + A n+1 At n+1 ) ■■■(! + A N+1 At N+1 )} 

I b N+1 [(1 - A N+1 At N+1 ) • • • (1 - ,4 n+1 At n+1 ) • • • (1 - AiAti)] ta 

-f^+^Adt'" 



The path integral for the isospin part thus becomes 

J dl% ■ ■ ■ dI?6(M N+1 I N+1 - I N )6(M N I N - /jv-i) • • • S(M 1 I 1 ~ h 



(4.14) 



oc 6(1 fTe 



■ft! Adt 



Ii). 



(4.15) 



We know that the solution of the classical equation I a + gfabcAtx^ = is 

I f = Te~ Jt * Adt h. (4.16) 

The unitary matrix U = Te~ Adt , which consists of real elements, has an inverse element 
[C/ _1 ] a fc = [CTfjofc = Uba- Thus, it can be seen that the solution of the parallel transport 
equation, Eq. ( |4.16| ) is equivalent to Eq. ( 4.15 ), that is, 



I f = Uh 

[u-Hff = [ii] a 

lb _ \T.m~-it! Adt -\ 



Uba[If] 



IfTe 



[h] a . 



(4.17) 



We see that the path integral agrees with the classical result. As a result it contributes 
as a constraint, the solution to the classical parallel transport equation. The path integral 
can be written in the form before the integration over the initial momentum Pf(ti) = j a 



dj 



f 

dX I ®h6(h-X) I ^xeK/*[ar* a -W e sJ( J / Te " Jii ^-'O 



(4.18) 



where A ac = gf a b c A b x^. Note that here we did not apply the gauge fixing for <j)Q from the 



subsection 4.1 because we will derive the Klein-Gordon constrained wave equation which 
has a gauge symmetry. In the next section and in Appendix |b|, the validity of this path 
integral will be supported by deriving the wave equation obtained as a constrained wave 
equation in the canonical quantization. In addition, Appendix [A| provides the derivation 
of the generalized Lorentz force equation Eq. ( p.6[ ) by varying the path integral. 
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5. Fundamental representation 



So far we have chosen the internal degrees of freedom of the classical particle to be in 
the adjoint representation. It is not difficult to write an action with the internal degrees 
in the fundamental representation. In the action of Eq. ( |2.3| ) we first write the structure 
constants f a \, c as —i(T b ) ac where T b are the generators. Then the 'charge vectors' I and 
J are taken in the vector space of the fundamental representation, so that we can replace 
fabcA b J a I c — > —iA b T b jJ l P. If the fundamental representation is real, for example if the 
gauge group is SO(N), we can choose the matrices T?- to be antisymmetric and purely 
imaginary. Then the Lagrangian becomes 

L = ^f~ \™ 2 + Jl ( p ~ igT«A«Pxn. (5.1) 

If the fundamental vector space is complex, for example when the gauge group is SU(N), 
the vectors / and J will be complex. Then the Lagrangian of Eq. ( |5.1| ) doubles the internal 
degrees of freedom. In order to avoid this doubling, we replace in this Lagrangian 

jaja _^ U i*p 

rf ahc A\P^ir*{-iT b 3 )A\l\ (5.2) 

where now the T a are Hermitian, but not all purely imaginary. Then we can write the 
Lagrangian as 



1 9 h 

— x 

2h 2 

We can then write the parallel transport equation for the complex charge vector i", 



L = -x 2 - -m 2 + il*\r - i<fl%A%P#>). (5.3) 



I* - igT^A^Px^ = 0, (5.4) 
and the generalization of the Lorentz force equation is now 

' " + ^(gAlT^P)-gT^P(d.A^ 



dt Vv^y <ft M ij 13 

= j t {h- 1 x,)+gK a F« u x» = 0, (5.5) 

where now we have defined K a = —Tf-P % P = K* a . The Lorentz force equation thus has 
the same form as Eq. ( 3.10| ), i.e. the same form for adjoint and fundamental representations. 
It can be easily seen that K a satisfies the same parallel transport equation as in Eq. (|2.7|), 

k a = -Tfitptp +r i P) 

= igP i P(T b T a - T a T b ) i:j A b ^ 

= -gf ahc A\K c x». (5.6) 

The constraints are similar to those in the adjoint representation, 

01 = P}, « 0, (5.7) 
4 = P\- iP l « 0, (5.8) 
4>h = Ph~ 0. (5.9) 
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We also note that 



P» = h~ L x^ - gAlK a . (5.10) 



After removing the variables I* 1 and Pj, , similarly to Eq. (4.7), the Lagrangian in the 
fundamental representation, Jzf^j, can be written down, 

= + p\i a + p} J* 1 + p h k-H 

= _ K n " H + ^ 2 " ^ m2 + p} ( p " igT ^ p±fi ) • (5 - n) 

And the path integral becomes 

J d\ J 9x< x 9P ll 9I i 9P}9h5[h- A][</>o,xoM0o]e^ d *^". (5.12) 

6. Derivation of the constrained wave equation 

In this section we derive the wave equation satisfied by the path integral. This links the 
classical formulation with the quantum one, in turn justifying the choice of the classical 
action. While we expect to find some generalization of the Klein-Gordon equation, the 
operator representation from the classical Hamiltonian brings ambiguity in ordering of 
operators [13]. We will see that the quantum Hamiltonian operator as derived from our path 
integral includes a term which does not remain invariant under the gauge transformations 
of Eq. (^^). We will verify this result by taking an inverse Weyl transformation of the 
classical Hamiltonian. The derivation is long, so for the sake of clarity we have gathered 
some of the intermediate calculations in Appendix |B|. 

The wave function at the final point (xt,If,tf) is the weighted sum of the wave func- 
tions at all possible starting positions (xj,/j,ij) , weighted by the kernel 

K(xf,If,tf;xi,Ii,U) = (qf,tf\qi,U} , (6.1) 

where K(xf,If,tf,Xi,Ii,ti) is the path integral Eq. ( |4. 18| ) and q stands for all variables, 
both internal and external. We are interested in the differential equation satisfied by the 
wave function. So we consider an infinitesimal evolution with £j = t and tf = t + At. We 
also write Xf = x, If = I and Xj = x — £, ij = I — rj. Then we can write the wave function 
as 

t + At) = N J *, J - 5 , / - t) K(x, t + Ai; x - 1 , / - „, «), (6.2) 

where N is a proportionality constant to be determined by matching wave functions in 
different times. Since At is infinitesimal, we can write 

(q f ,t + At\q h t) = (q f \e-T> HAt \ qi ) 



J dp e i (9/-^)p/^ e -|^(P.(9/+^)/2)At 

e** M , (6.3) 
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where H is a Hamiltonian operator, H c (p,q) is a Hamiltonian function and is a La- 
grangian in a discrete time interval. Here p stands for the momenta canonically conjugate 



to all variables. We identify the Lagrangian Jzf from Eq. ( [1.18 ) and expand Te 
first order in e = AAi, using [A] ac = gfabcA^x^ to write 



to 



^At 



2A 



(Xf - Xi)' 



m 



2e 



I b f {5 ba -ghcaAl{x f -x t Y)-It +0(e 



(6.4) 



where r) = If — Ii and £ = (x/ — X{) . According to the Weyl correspondence [13] the 
classical Hamiltonian H c has to be evaluated at the midpoint 

(q j +l\e- l ^ +1 ~ t ^ /h \q j } oc (q j+1 \ [ e- iit ^- t ^ (p ' q)/h A{p,q)dpdq\q j ) 



oc / dpe 



(*3 + l— *j) 



(6.5) 



with 



A(p, q) = j due iq - u ' h \p - u/2) (p + u/2\. (6.6) 
Thus the isospin dependent part of the action takes the midpoint value, denoted by a 'bar', 

Si = P?nK ~ In-l + gfa b XJi{x» - X^)}. (6.7) 



Using the expression for the path integral from Eq. ([4.18] ) in Eq. (|6.2| ) and using Eq. (|6.4|), 
we can write the wave function at t + At as 



with 



where 



i/>(x, t + At) = N J dr] J d&(x -£,I-v, t)K{£, V , t), 

J-oo JO 



(6.8) 
(6.9) 



1 



A* x = x 1 + M* - x (i - vW 



A^x) x / - -^d v A„(x) x / - -A^x) x v 



e + o^ 2 7])- (6.io) 



Up to the second order in £ and the first order in rj, the exponent in Eq. ( |6.9| ) can be 
expressed as 



i 

h 

i ( 1 



|- +f u ll xi--A li xT,\ e- ^f{d^A v x irec + jv - ^ T 



(6.11) 
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where we have defined 



a ^ = G ^f[ Av xI--A v xrj\ a , 

e = e + «a 



(6.12) 



and G^ u is the matrix inverse of G^ , so that G^G uX = 5* . On the right hand side of 
Eq. (|6 . 2|) , we expand the wave function ifj(x — £, I — 77, t) up to the second order in £ and 77, 

- e, / - T7, = /, *) - - + e 

+ ^f^W + ^V^d^. (6.13) 

Inclusion of higher orders £ and 77 in a wave function expansion brings higher order contri- 
bution in e, as we will see below. Now we are ready to integrate to get the wave function 
at t + At, 

il>(x, t + At) = N J dXdjd^di [i> - (e - ea^)d^ - rfd^ + - ea») V a d a d^ 



xe'2£ 



{j-J^G^C ■-iiaPG liv a v -ie?%-+ij a r) a ]/h 



(6.14) 



Let us first perform the Gaussian integral for £. Odd order terms in £ vanish, so those are 
dropped and only the terms up to first order in e are kept in the result, which reads 



V>(M + At) = TV / dXdjdT]d^I^e [ - l 2 



(6.15) 



h is defined as 



x e 



1 



In this equation we have written 

Ar-i (i ) = J d&irJ 11 ^ = (27^)^/ 2 (detG^)- 1 /2, 
J di^i v e^ G ^" =iheN- 1 {j)G» u . 



(6.16) 



.17) 
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As mentioned earlier it can be seen that higher orders in £ bring higher orders in e. A 
similar statement is true about rj. Next for integration for r] and j let us keep terms only 
in the first order in e, 



(j • (A, x I)) 2 + -(j ■ (A, x v )) 2 - j ■ (A, x I)j ■ (A" x r,) 



where we have used the fact that G^ u = rj^ + 0(e). Let us write the result as 



.18) 



ip(x,t + At) = N J dXIe~ le ^ 
o 



where 



j djd V Ii:e- L ni a " G ^ + L n^ 
djdvN-^j) 



eh 



./■>i 



+e / djdrjN-^j) [-a^d^ + a" V a d a d^} 



(j ■ (A, x I)f + -(j ■ (A, x V )f 



-e— j djd V N-\j) 



-j ■ (A M x I)j ■ (A" x v )} 



1 



J-'l 



ihe 



djdriN-'^G^dMen^ 



h + h + h + h- 



(6.19) 



(6.20) 



The integrals - ■ ■ , I4 are calculated in Appendix |B|. Using the results given there, we 
can write the final result 



/ = AT-i(O)^ + AiV- x (0)^ [-(dpA" x I) a d a ^ - 2(A" x I) a d a d^ 

+\tt{A ■ A)^ + (A" x I) b facbA c ^ + (A„ x I) a (A» x I) b dM + d 2 ^ 



,(6.21) 



where we have brought back A using e = A At, and written Tr(A • A) = j acb^hdaA^ . If 

_■ «£ 

all the terms, including the mass term e le 2 ; are now gathered, we find the wave equation 
with respect to the parameter t. Let us first choose the proportionality constant N such 
that 

[CO [CO 

1 = N dAiV _1 (0) = N(i2hAtn) D/2 / d\X D/2 
Jo Jo 

(i2hAt7r) D / 2 , n/9 , , 

= lim N- ; - A D / 2+1 , 6.22) 

A^oo D 2 + 1 ' v ; 
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where A is to be sent to infinity. That is, 

D/2 + 1 1 



N 



{i2hAtTT) D / 2 A D / 2 +! ' 



(6.23) 



Going back to Eq. (6.21), we can write the wave function t/j(x, I,t + At) at t + At as 



ip(x,I,t + At) - ip(x,I,t) 



-(SLA" x J) a d */; - 2(A» x I) a d a d^ + -Tr(A • A)^ 



to 



+(A» x IfUaA'^i, + (A„ x I)'(A" x /)°a„W + - -p-</> 



2 Jo 



x I) a ^-PjV - 2(A" x /) a Q ) P/P^ 



1 



+ -Tr(A ■ A)1> + (A» x I)"f acb Al-P^ 



n 



, (6.24) 



where we have used that d 2 = -<9 2 + d 2 = -p-(-P 2 ) - ^P 2 = — ^ 
The wave equation with respect to the worldline parameter t is 



1 P 2 and P? = %di*. 



ih 



where 



-N 



| /'OO 

--N 
2 ./q 



dAAiV -1 ^) -ih(d^ x /) a P/ + 2(A^ x I) a PfP^ 

+^Tr(A • A) + x I) b f acb A^Pf 

— (A^ x /) a (^ x l) h FfP$ - P 2 -to 2 

dAAiV -1 (0) [-(P 2 + to 2 ) - 2{Alk a )P^ 



» v (A a A b v )k a K b - ih^Af 1 x I) a P 7 a 



-77 



/i 2 . 



+^Tr(A • A) + ih(A» x I) b f acb A^Pf 



1 /-co 

-JV / dAAA^ 1 (0)[L>V> + A^], 

2 Jo 



(6.25) 



6 rc 



#° = -AacPf/' 

£ty = [(P„ + ^A^) 2 + m 2 ]QQ...ppV, 



iT^A" x I) a Pf - %Tr(A ■ A) - ih{A» x I) h f acb A c Pf 



(6.26) 
(6.27) 

(6.28) 



We will see below that the operator D + A is an inverse Weyl transform (the Wigner trans- 
form) of the classical Hamiltonian H = |[(P M + gA^K a ) 2 + to 2 ]. The notation QQ ■ ■ ■ PP 
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means that all the momentum operators are put on the right while the position operators 
are on the left. It can be seen that the right hand side in Eq. ( |6.25| ) is divergent, since 



N 



roo poo 

/ dXXN-^O) = N{i2hAt^) D/2 / d\X D/2+1 
Jo Jo 



D/2 + 1 1 (i2hAtn) D / 2 D/2+2 
(i2HAtTt) D / 2 A D / 2 +! D/2 + 2 
D/2 + 1 



D/2 + 2 



A. (6.29) 



Eq. (|6.25|) can be written as 



ih^(x»,I a ,t) = I^Z|±1a(DV + A# (6-30) 

Dividing the both sides by the divergent factor and taking the limit A — > oo, we find 

DiP + A^= lim ih D Jl + 2 \?-Mx^I a ,t) = Q. (6.31) 
A-s>oo D/2 + 1 A at 

This is the wave equation for tp. For a fundamental representation, the wave equation is 
obtained by replacing, e.g., Pf f a b c A b ^I c by —iP\Tf-A b ^P. Note that in this substitution 
one has to maintain the order of indices in /a& c . Then in the fundamental representation 
Eqs. fls!26D- (ET28D are replaced by 



K a = iT^PjP, (6.32) 

(6.33) 

if). (6.34) 



-ij 

_L nr, A ar T a TO pi\ 2 j_ ^2, 

\QQ-PP 

r , .... ft 2 

Aip 



D^ = [(P„ + igApfjPP}) 2 + m%A P p^, (6.33) 



hg{d,A b ^P)P} + ^(T^A^) + ihg 2 {A^T«P)T b kl A b ^ 



6.1 Operator ordering 



In this subsection we verify the wave equation of Eq. ( 6.31| ) by using a general mathematical 



formula which relates the quantum Hamiltonian operator to the Hamiltonian function in 
the path integral. The map from the operator to a function in the phase space is called 
the Weyl transformation. The inverse map is called the Wigner transformation. The 
correspondence between the two spaces is one-to-one. The Wigner transformation from a 
function a(p, q) to an operator A(P, Q) is given in a compact form as [13, 14]. 



A(P,Q) = [e"*>-*a(p,q)]y_ ¥ p t9 _^. QQ ... rr . (6.35) 

where QQ ■ ■ ■ PP means that all the momentum operators are placed on the right side 
while the position operators on the left side. 

The classical Hamiltonian H used in the path integral was ^[(P /Ji +gA'/ l K a ) 2 +m 2 }. The 
Wigner transformation of that function on the phase space should give us the appropriate 
operator to be used in the Schrodinger equation. We note that the terms in H for which 



- 18 - 



there is an ordering ambiguity are 2P^ l A^ 1 K a and rf v A^K a A b K b . Plugging these functions 
into Eq. ( 6.35 ), we get 



h d a 

e 2i aPfj. dxv (2pv- A a K a ) 



2AKP + -dpAUK 



(6.36) 



and 



9 ° a (^ u A*K a A b u K b ) 
(An x I) a (A^ x I) b PfP b j 

' Ab ^^ {facdfbefPfldP!lf) 

_l. a A b^_d_d_ 9 d_ 

2 ' 4 dPf dia dPf dl h 

(A„ x I) a (A^ x ifPfP} 



(facdfbefPfI d PfI f ) 



±{-2A^f bca (A» x lf)Pf - ~2(f acb A c ) ■ (f bda A d )) 

h 2 . 



(6.37) 



(Afj, x I) a (A» x I) b PfP b - ~A c J bca {A» x I) b Pf - '-TriA ■ A). 



b £>a 



(6.38) 



The additional pieces, ^d^K, -f^/ 6ca (^ x I) b Pf and -^Tr[A • A], exactly match 
with the terms in Eq. ( 6.28 ), the expression for Aijj. Let us rewrite the expressions ( 6.36j ) 
and ( |6.38| ) in a different form. The operator ( |6.36D can be written as 



,» ep M afp 2PAK (P • A a + A a • P)K a . 



(6.39) 



Also, using the relations 



\{A,xI) b ,Pf] = gf bcd A c fl [I d ,P?] 



(6.40) 



and 



x I) a (A» x ifP/P* = (A„ x I) a [(A^ x /) 6 ,P/]P / 6 + (A„ x /)°P/(^ x I) 6 ^ 
= (A„ x I) a (ihf bca Al)P} + x /) a P/(^ x /) 6 p) 
= (A„ x 7) a (^/ 6ca ApP 7 6 + (^rA^'i 6 ) 
= (A„ x I) b (ihf acb A^)Pf + (^A-i^*) 
= -(A„ x I) b (ihf bca A^)P? + (^iTM 6 ^ 6 ), (6.41) 

we can rewrite Eq. (6.38) as 



, 2 R i a? a (jf v A a K a A b v K b ) 



■(A„ x I) b (ihf bca Al)Pf + (A a k a )(A b »K b ) 



+ihA c fl f bca (A» x /) ft P/ - T Tr(yl • A) 



{Alk a A bll k b ) - — Tr(A ■ A). 



(6.42) 



- 19 - 



Thus the Hamiltonian operator can be written without assigning a particular ordering, 



H = h - 
2 



(P„ + gA a K a f + m z - g 2 '-Tr(A ■ A) 



h 2 



(6.43) 



We clearly see that the Hamiltonian operator consists of the operator naively replaced by 
the classical Hamiltonian and the gauge non-invariant term g 2 ^Tr(A - A). As far as the 
operators orderings are relevant, there is no reason to expect that the Hamiltonian function 
in the path integral should correspond to operator in which the phase space variables have 
been naively substituted by the corresponding quantum operators. We note that the last 
term has been previously found in the literature [12]. 

The Hamiltonian commutes with the internal angular momentum K 2 , so one can factor 
out the eigenfunction for K 2 . Because of using real and bosonic variables, K a behaves as 
an ordinary angular momentum and hence has integer eigenvalues. 



7. Conclusion 

In this paper we have considered the question: What is a classical non-Abelian point par- 
ticle? More specifically, what is the classical dynamics of such a particle in a background 
non-Abelian gauge field? We started from a classical action describing the position of the 
particle as well as its charge, described by a dynamical 'internal' vector in some representa- 
tion of the gauge group. We found that when this internal vector is in either the adjoint or 
the fundamental representation, the charge vector that enters the generalized Lorentz force 
equation is an adjoint vector constructed from the original charge vector and its conjugate 
momentum. So the charge vector that determines the space-time trajectory of the particle 
is in the adjoint representation in both cases. 

The equations were originally derived from quantum theory, so we decided to quantize 
the classical action in the path integral formalism as a kind of cross check. Using the 
worldline formalism so as to include particles of zero mass, we found that the sum over paths 
includes only those paths along which the internal charge vector is parallel transported. 
We also derived a wave equation from this path integral and showed that the Hamiltonian 
operator in the wave equation exactly matches with the Hamiltonian operator transformed 
from the classical Hamiltonian function by the Weyl correspondence. 

There are however some differences between the quantum theory of a non-Abelian 
charged field and that of a non-Abelian point particle as constructed from the classical ac- 
tion. That is not a failure of quantization, nor a contradiction with established knowledge. 
The first point of departure is the fact that the classical 'isospin', the charge vector in the 
action, is a continuous variable, whereas for quantum particles, the isospin / is quantized, 
with fixed I 2 . This can be resolved in the following way. The charge vector K in our 
setup corresponds to the isospin vector operator upon quantization. By construction the 
vector K is like an angular momentum operator in the internal vector space, so it will have 
discrete eigenvalues when the theory is quantized. The quantum Hamiltonian commutes 
with K 2 as can be easily checked, so a particle in a given isospin eigenstate remains in that 
eigenstate. Then the quantum particle with a fixed isospin corresponds to the particle in an 
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eigenstate of K 2 . However, since K corresponds to x x p in the internal space, it can have 
only integer eigenvalues. This is analogous to there being no truly classical description of 
half-integer spin. The solution would be to use anti-commuting variables for the internal 
charge vector. That would allow, for example, half-integer isospins when the gauge group 
is SU(2). 

Another discrepancy is the appearance of the term ^<? 2 7i 2 Tr(A • A) in the quantum 
Hamiltonian operator, Eq. ( |6.43 ). This term breaks the gauge symmetry present in the 



classical Hamiltonian, although the Hamiltonian is still symmetric under constant internal 
rotations. This term appears to be a genuine effect of quantization: the Hamiltonian 
operator as derived from the path integral exactly matches with the Hamiltonian operator 
constructed from the classical Hamiltonian function by the Weyl correspondence. While 
a Hamiltonian constructed directly from a non-Abelian gauge symmetric quantum field 
theory is not expected to contain such a term, we note that the anomalous term is not 
unknown in the literature [12]. 

The source of this discrepancy is the following. The constraint which implements gauge 
transformations comes from the Hamiltonian for the gauge field, but we have treated the 
non-Abelian gauge field A^ as a background field, ignoring its dynamics. We could of 
course try to include the Lagrangian for the gauge field. But it is known that for an 
electrically charged point particle in a background electromagnetic field the joint action 
leads to inconsistencies, stemming from the fact that the field due to the charged particle 
itself diverges at the position of the latter [15]. The problem is resolved by starting from 
the Lorentz-Dirac equation instead of the Lorentz force equation, and making further 
modifications so that quantization leads to the Dirac equation [16, 17]. 

The appearance of the gauge symmetry breaking term in the quantum Hamiltonian for 
the non-Abelian point particle is related to this classical difficulty of defining the action of a 
point particle in a dynamical gauge field. In the Hamiltonian picture, the constraint which 
implements gauge symmetry appears only if the gauge fields have their own dynamics. 
But if the non-Abelian point particle is coupled to a dynamical gauge field, the radiation 
reaction must be included as for the ordinary electric charge. In that case, Wong's gener- 
alization of the Lorentz force equation will have to be replaced by something analogous to 
the Lorentz-Dirac equation, and a corresponding action, as the starting point. 
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Appendices 

A. Derivation of the equation of motion 



In this appendix we show that by varying the path integral in Eq. ( [4.18 ) with respect to 



X 1 *, we can recover the classical equation of motion for x^ from the vanishing first order 
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variation, 

t f 

/•OO rOO p 2 — f Adt 

= 6 dj d\@h5{h-\) Stx^efJ^^-^MiVfTe K -A), 



(A.l) 



J-oo JO J 

where A ac = gfg^Ax^. Using the relation 

'to 

where U(tb,t a ) = Te~f*a Adt , we can write 



SU(t b , to) = - [ " dt'U{t b , t')5A(t')U(t', t a ), (A.2) 



d -U(t f ,U) = - [ 1 dt'Uitf^-^—A^Uit'.U 



' dt'U(t f ,t')[S(t' - t)d„A u {t')x v (t') - A^S(t' - t)]U(t',U) 
-U(t f ,t)drA v (t)d?U(t,t a ) - j t [-U(t f ,t)A^U(t,t i )} 

-u(t f ,t)[dpMt)#W(t,ti) 

+U(t f ,t)A u ± u A fM U(t,t i ) + U(t f ,t)d u A^x v U{t,ti) 
-U(t f ,t)A f ,A u x I/ U(t,t i ) 

-U(t f ,t)F^ u ± ,/ U(t,t i ), (A.3) 



where [F^ u ] ac = gfabc[d,iA b v - dyA^ + gf bde A^Al\. Recalling that it has been defined in 
Eq. ( 4.18 ) that j a = Pf{ti), the relevant term jlf g x ji/ t \ U(tf, ti) can be written 

P?(ti)[-I f U(t f ,t)F^U(t,ti)] a = Pfit^-I^F^Uit^a 

= [-I(t)F^x v U(t,ti)] a P?{ti) 
= -I(t)F^± u Pj(t) 

'bc& fit/ 



-I a {t)gf abc F%x»Pf{t) 

-gK b F*V, (A.4) 



where K a = f abc PjI c . Therefore, the vanishing first order in 5x^ in Eq. ( |A.1| ) leads to the 
equation of motion for x^. 

J dt{-j t (h'%) + J? [-IfU{t f ,t)F^U{t,U)]aS^ 

+ 9 K"fL x" = 0. (A.5) 



d ( mx M ^ | ^„ 



dt VV^F 



(A.6) 
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B. Derivation of the constrained wave equation 

The kernel was written from Eq. ( |6.11| ) in the following notation, 

where 



oo poo 

dj / dX exp i - 

-oo JO I n 



1 — — f Til 

-eG^--a»G^ + jr ] -e— 



(B.l) 



= r,^ - -ef[(d^A v x If + (d v A„ x I)% 
a» = G^j a [A u xI-^A u x v } a , 

e = e + 



(B.2) 



The wave function ip(x, t + At) at t + At is a sum of the wave functions ip(x — £, J — 77, t) 
with the weight K(£,7],t). 

i/j(x, t + At) = N J dij J d£i(,(x -CJ-V, t)K(£, rj, t) 

= n J dxdjdr/dt U - o^e - 9 a ^v a + \d^ee 



+ -d a d b tf;rj a V b + <9 a d M V£V 



N f dXdjd^dl [if, - d^(e - eaP) - d a ifnf + d a 3^{e - ea^rf 



iV / dXdjdrjd^I^e 



x e2R. 



(B.3) 



which was Eq. ( 6.15| ). Let us perform the Gaussian integral 1^ first. In the end we are only 
interested in the first orders in e after suitable choice of normalization factor N. 



+ -dad b ^ a V b - 8 a d^(e ~ ea")rj a 



1 



e 2hc 



N-Hj) U - d a ^rf + -d a d b ^ a rf> ) + N-\j)e[d^{a> 1 ) - d a d^{a»)r, a ] 



(B.4) 



- 23 - 



where it has been defined 



N'^j) = J dle^ G ^ v = (^evr^detG^r 1 / 2 . (B.5) 



The vanishing exponential integrals with odd multiples in £ have been dropped out and us- 
ing the formulae <5(det G^) = (det G^G^SG^ and bG^ v = -G^G V0 ' bG pa , the following 
integral has been evaluated. 



/ 



dCeCe^ G ^ = -i2he 



dN-^j) 



dG 



-i2he(i2heir) D/2 



dG- 1 ' 2 
dG 



[IV 



-i(detG^)- 3 / 2 



= -i2he(i2heTr) D ' 2 
= iheN' l {j)G^. 
Next what remains is to integrate over (rj,j). 

1 = J djd^e-^G^H^ 



(det G^GV 



(B.6) 



(B.7) 



Since only the first order in e is needed, in the exponential e 2ft aMG V a " on iy the terms up 
to first order in e are kept. 



1- —j [A^ x I- -Ap x v ) CTj A„ x I- -A v x v + 



l-—j[A li xI--A li Xf]\ rf v j \A V x / - -A v x r, + 0(e 2 ) 









2h 




2fr 








+0(e 2 ), 



e r 



.7'/ 



1 

4- 



B 



where G^ v = r/^ + O(e) has been used. The function e^ ri is proven to be useful since a 
function of rj can be converted to a differential operator with respect to j. 



J djdT]e^ T >f(r])g(j) = J djdrjg(j)f 

-I 



h d 
i dj 



g(j) 

= (2nh) n J dj6(j)f (-^ 9 {j) 
= (^) n f{~§ ] )9(j)\ J =o, 



(B.9) 



where n in (2irh) n is dimension of j a . Note that we will omit (2irh) n in later calculations 
since it can be absorbed into the normalization factor N . It is worthwhile to recognize two 
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properties for calculational convenience. The first is that to get non-vanishing terms the 
number of powers in rj must be the same as that of j's. If k ^ m, 



djdr]ei jv r] k j m = 



(B.10) 



The second is that nth derivative of ./V (J) with respect j brings e n order. For instance, 



dN~ l (j) 



(i2eh<K) D / 2 
1 



i(det G^ 2 



Qja 



-N^ii^CT-^- = -N- 1 (j)G> iU B? l 
2 dj a 2 



(B.ll) 



It can be easily seen that the further derivatives also have the similar properties. This 
implies that higher than first derivatives of iV _1 (j) with respect to j can be discarded 
because we want to keep terms only up to the first order in e. The final integral I is 
written in four pieces. 



1 = J djd^e-^^+Ti™ =I x +I 2 +h + I 4 , 



(B.12) 



where 



h 
h 
h 



j djdriN- l (j) U - d a ^r? + i&^y) e^", (B.13) 
e J djdqN-^j^ian ~ dad^a^ei^, (B.14) 
~ J djd^N-'U) (j {A^ x l)j (A» x /) + ~j x v )j (A» x r,) 

(B.15) 
(B.16) 



X ( tj} - d a ^ a + -d a d h ^ a rf ) e~^\ 



Let us evaluate the first integral I\. 



1 



djdvN-'ij) $ - doM>rf + -dad^rf ei« 



h d 



djdrjN (j) ( ij) - d a ip-— ] en 



(B.17) 



where ^d a dbipr] ci "rj b has been dropped out since it corresponds to the second derivative of 
A^ _1 (j) with respect to j, which is in the second order in e. Thus, using Eq. ( p.ll| ) 

h = N-\<d)i, + d a ^-^- a N-\j)\ j=0 
i oj a 



N-^Qty + da^N-^CTB^ 
ihe 



(d^ x i) a d a ^ 



(B.18) 
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Next, in I2 we discard jrj flu (A I/ x /) and ^JV^i^u x v)v a since they do not have the same 
powers in j and 77, so those terms vanish or contribute to higher order in e. 

h = e [ djd V N-\j)[dM) - dadMa^jei™ 



dMGT [A u xl--A u xri 



e J djdrjN-^j) 



dad^jG^ [ A v x / - ~A„ x 77 ) rf 



e'> 



.IV 



e f djdrjN- 1 ^ 
+0(e 2 ), 



.1 



( -j-A u xr/)- d a d^(j (A, x I))rj a 



eh 



.in 



(B.19) 



where G^ has been replaced by rf 111 to keep only the first order in e. And the first term 
in Eq. (|B.19 ) does not contribute, i.e., 



djdrjN^iJWd^ (-j Qa, xr,yj el* 



-e J djdrj 
+ O(e 2 ). 



^UyrTd^— Qa, x /f + Qa„ x -^(N-Hj)^) 



r.l'l 



(B.20) 



Finally, the only non- vanishing part in I2 is 



h = ~ej djd-nN-^dednWV {A v x l))rfe^ 
= e J ^driN-^dad^ — U (A v x I))ei 
= e f djd^ijW^-iAv x Ifei™ 



eN- l (Q)d a d^-{A^ x If. 



(B.21) 



By the same reason as for I2, in ^3 only for the non-zero first order term in e, only terms 
in the same power of 77 and j apart from N~ l (j) are kept, 



2h 



J djdrjN- 1 (j) (A„ x /) j (A» x /) + i j (A„ x 77) j (A" x 77) - j (A^ x /) j (A" x 77) 



x l v - 9»W + -d a d b ^ a v h ) e**> 



^31 + -^32 + ^33- 



-j (A„ x 77) j (A» xrj)i)-j (A^ x /) j [A 11 x 77) (-d a 4>r] a ) 



+j (Ap x J) j (A" x /) ( -dad^r? 



en 



j' 1 



(B.22) 
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We calculate J3 in three parts, /31, I32 and I33. 

ie f 1 i 

hi = -2hJ djdvN-^j)-^ (A^ x n)j (A* x n) e*" 

= ^N-\U)\{f acb A^t V {heaAl)i> 

= l -^-N- 1 (O)Tr (A ■ A)if). (B.23) 
8 



7- - K 



^ I djdr / iV- 1 (j)[(-j (A„ x (A" x 7/))(- W)]e*™ 

f djdr]N-\j)(-d a ^e^^^ b [(-j (A,x I))(-A^ x j) b ] 



= ~2ft _ 

= ^iV- 1 (0)(-9^)[-(^ x I)V u (+fac b Al)} 

= l -^N-\0)(A» x I) b f acb A^d a i>. (B.24) 

= ~ / ^iv- 1 ^) Q w) 7^7^ x (A" x /))] 

= -^AT-^O) QoAv) x J)°(A" x /) 6 

= ^AT-^O)^ x I) a {A» x /) b a a ^V. (B.25) 
The last integral I4 is straightforward. 

h = ! Tj tiW^G^d^eV" = % ^N-\0)^. (B.26) 
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